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SYMBOLS 


ai rate of speed change in the first segment of a three-segment speed 

profile 

a 3 rate of speed change in the third segment of a three-segment speed 

profile 

Aa in speed profile, acceleration, miraerically positive 

in speed profile, deceleration, numerically negative 

C initial circle of imit radius 

o 

final circle of unit radius 

d^ signed distance of Pq^v’ positive if the corresponding vector ri^ 

points away from 

d^ signed distance of Pv^f 9 positive if the corresponding vector n^ 

points away from P^ 

d signed distance from P^ to another point P which has the same 

^ subscripts as d 

g acceleration of gravity 

h^ cruising altitude 

h^ final altitude 

h rate of altitude change; the sink rate if negative 

k numerical constant for path stretching, 0 < k < 1 

I normalized total path length of horizontal tra^jectory 

5,^ normalized total path length between and in problem (d) 

L as an algebraic quantity, the total horizontal path length; as a 

maneuver symbol, left turn 

Li minimum distance to be traversed in time interval t when the speed 

is constrained to lie between and 

L 2 maximum distance to be traversed in time interval t when the speed 

is constrained to lie between and 

desired path length for stretching 


V 


I 
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Lg stretched path length 


L . 
min 


L 

max 


n^ 


n 

P 


minimum distance to he traversed in time interval t when the speed 

is constrained to lie between V , and V 

min max 

maximum distance to be traversed in time interval t when the speed 

is constrained to lie between V . and V 

min max 

unit vector having the same direction as the initial heading 
unit vector having the same direction as the final heading 
initial horizontal position of the aircraft 


final horizontal position of the aircraft 

P^ point of intersection (vertex) of the two lines containing the two 

unit vectors and n^, respectively 


P switching point; first subscript is either o or f, referring to 

^ points in the initial line and in the final line , respectively; the 

second subscript distinguishes different switching points in 
problem (a) 


R as an algebraic quantity, turning radius; as a maneuver symbol, right 

turn 


R . minimum turning radius 

S as a maneuver symbol, straight flight 

i = 1, 6; in problem (b), switching lines; defined in equations (30) 
through ( 35 ) and figure 12 

t time; in the controlled time^of-arrival problem, time to go before 

reaching the final point P^ 

ti in the speed profile, time when constant-speed segment begins 

t 2 in the speed profile, time when constant-speed segment ends 

t . in the controlled time-of-arrival problem, earliest arrival time 

min 

achievable with speed control alone 

t^ax controlled time-of-arrival problem, latest arrival time 

achievable with speed control alone 

At 2 in the speed profile, duration of time when the speed remains constant 

Atj^ in the altitude profile, duration of time when the aircraft is 

descending 
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V 

Vo 

V 

max 

V . 

min 


V 

up 

X 

X* 


y 

y* 

^max 


ground speed 
initial ground speed 
final ground speed 

lowest achievable speed within the constraints of t, 

maximum ground speed of the aircraft in the terminal area 
minimum ground speed of the aircraft in the tenninal area 
constant speed in the three-segment speed profile 

highest achievable speed within the constraints of tj Ag^, A^ 

abscissa of Pq in a coordinate system where the final point is 

taken as the origin 5 the direction of n^^ as the positive x-axis , 
and the direction to the right of n^ as the positive y-axis 

normalized x, where the minimiim turning radius taken as the 

lonit of distance 

in problem (d) , the distance of the aircraft along the first segment, 
measuring from the point of intersection 

in problem (d), the distance from the point of intersection to a 
desired final point in the second segment 

ordinate of Pq , measured in the same coordinate system as x 

normalized ordinate of y, where the minimum turning radius ^\nin 
taken as the unit of distance 

bank angle 

maximum bank angle 

heading of aircraft 

initial heading of aircraft 

final heading of aircraft 

total heading change in the initial turn 

heading of the straight-line segment of the horizontal trajectory 
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TERMINAL-AREA GUIDANCE ALGORITHMS FOR 


AUTOMATED AIR-TRAFFIC CONTROL 
Heinz Erzberger and Homer Q. Lee 
Ames Research Center 


SUMMARY 


In present air-traffic control systems , certain terminal-area guidance 
problems are solved manually by the pilot and the controller. The solution of 
these problems could be automated in a new generation system with the use of 
ground and airborne computers and an improved navigation system (such as the 
proposed microwave landing system). The objectives of this study are to for- 
mulate these problems analytically, and to obtain solutions to them in the 
form of computer- oriented algorithms. 

The principal requirement is to construct a flyable , three-dimensional 
trajectory that begins at the current aircraft position, heading, speed, and 
altitude, and that terminates at a prescribed position, heading, speed, alti- 
tude, and time. The terminal position is normally the Instrument Landing Sys- 
tem (ILS) gate, or a waypoint, and the terminal time is the assigned landing 
slot. Such precise time-position control, "four-dimensional guidance," is 
proposed as an advanced air-traffic control technique. The algorithms devel- 
oped in this report are applicable to all possible combinations of initial and 
final conditions, and thus can be us^d in a closed-loop feedback law. 

The construction of the required trajectory is divided into three parts. 
First, the horizontal trajectory is calculated with a constraint on the turn- 
ing radius. Algorithms for constructing the horizontal trajectory are given 
for each of four criteria: (a) to fly to a waypoint and arrive there with a 

specified heading, (b) to fly to a directed line in space such as a VHF Omni- 
Range (VOR) radial, (c) to fly to a waypoint without a constraint on the head- 
ing at the waypoint, and (d) to change smoothly from flight along one line or 
VOR radial to another that intersects it. For each criterion the trajectories 
are synthesized from straight lines and portions of circles. Assignment of 
the trajectories is based on minimizing the length of the path so as to reduce 
maneuver time and conserve airspace and is given in the form of switching 
diagrams . 

Second, the speed profile is calculated based on the length of the hori- 
zontal trajectory and on the desired time over the waypoint or ILS gate. Con- 
straints are the minimum and maximum airspeed, the acceleration and decelera- 
tion of the aircraft, and the desired airspeed at the terminal time. The 
algorithm first determines if a speed profile satisfying these constraints 
exists. If it does, the speed profile is synthesized; if it does not, a tech- 
nique is used to stretch the path length of the horizontal trajectory by the 
amoiint required to synthesize a flyable speed profile. 



Third, the altitude profile is synthesized so that the aircraft maintains 
the approach altitude as long as possible before descending to the specified 
final altitude at the waypoint. Input to this calculation is the length of 
the horizontal path, the desired descent rate, and the length of the final 
deceleration interval. The descent to the specified final altitude is timed 
so it does not coincide with the deceleration interval. 

Finally, the completely specified flight profile is arranged into a time 
sequence of guidance commands that can be used as inputs to a flight director 
or autopilot. 


INTRODUCTION 


The design of aircraft navigation and flight-control systems has long 
followed a trend toward automation of the numerous guidance, control, and nav- 
igation functions once performed manually by the flight crew. Recently, this 
trend has accelerated and expanded in scope, with the momentum behind it being 
provided by the development of flight-worthy computers and by the increasing 
complexity of modern aircraft and of the air-traffic environment. Complex 
problems in all phases of aircraft operation indicate that such automation is 
not a luxury, but rather a prerequisite for safe and economic flight. Fortu- 
nately, the timely development of airborne computers has provided the means to 
achieve the degree of automation required to alleviate these problems . 

Concepts have been proposed for flight-control systems capable of flying 
an aircraft automatically from takeoff to landing (refs. 1 and 2), with the 
pilot purely in the role of systems manager. Although such completely auto- 
matic systems are now under study, this study focuses on the problem of auto- 
mating a segment of the flight during which pilot workload is particularly 
heavy; the segment is within the terminal area, roughly within a 50-mile 
radius of an airport. For a subsonic jet transport on a landing approach, it 
corresponds approximately to the final 15 minutes of flight. 

Pilot workload in the terminal area is heavy as a result of deficiencies 
in the existing flight-control and air-traffic control systems. In this area, 
the aircraft often must be flown manually, even though an automatic landing 
system and an autopilot are on board. The automatic landing system can be 
used to fly the aircraft only along the final landing approach, after the air- 
craft is at the proper altitude, position, and heading to receive the Instru- 
ment Landing System (ILS) glideslope and localizer signals. Similarly, the 
standard autopilot has limited use in the earlier portions of terminal-area 
flight because it was never designed to provide precision flight-path control 
along the curved, decelerating, and descending flight path which is frequently 
reqiiired before the final approach. 

The guidance of aircraft in busy terminal areas, especially in weather 
conditions requiring instrument flight , is conducted primarily by air-traffic 
control. Controllers monitor the position and movement of aircraft on radar 
and issue guidance instructions to the pilots, who close the guidance loop by 
executing them as quickly and faithfully as possible. For landings, the 
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controllers direct the randomly arriving aircraft into the final approach to 
achieve a high landing rate while maintaining a safe separation distance. The 
task is complicated by the differing airspeeds (and other performance charac- 
teristics) between aircraft on a common path. 

The pilot and controller workload of this system of terminal-area guid- 
ance has been studied for a number of high traffic terminal areas (ref. 3). 

For arrivals into JFK International Airport, pilot-ATC commimication requires 
as much as 50 percent of the last 10 minutes before touchdown. During this 
time, the pilot performs an average of l8 communications and navigation opera- 
tions, a workload which pilots consider excessive. 

To reduce the pilot-ATC workload in the terminal area, an automatic 
flight-control system must provide precision flight-path control along the 
specified trajectory. In addition, the system must also maintain synchroniza- 
tion of the aircraft position with a specified time schedule; this requirement 
is especially important for efficient air-traffic control. All methods pre- 
viously studied particularly stress the importance of maintaining accurate 
time-posipion control of aircraft in the terminal area (refs. H-6). Accurate 
time-position control of the aircraft should reduce landing delays due to 
holding and increase landing rates; however, it cannot be achieved with exist- 
ing air-traffic control methods and equipment. 

What techniques are available for guiding an aircraft automatically along 
a specified time-position trajectory? The technique suggested in reference 1 
is to precompute the trajectory in a ground facility, store it on tape, and 
enter it into the on-board system before takeoff. The on-board equipment 
would consist of a means for storing the precomputed trajectory and a modified 
autopilot for flying the aircraft along it. Although this technique requires 
only relatively simple on-board equipment, it suffers from severe practical 
disadvantages - the dependence on a ground facility and the difficulty of 
making fast trajectory changes in flight. For some segments of a flight, such 
as climb to and descent from cruise of a supersonic transport, these disad- 
vantages may not be significant. In terminal-area operations, however, even 
with precise scheduling of aircraft, unexpected trajectory changes can never 
be completely eliminated, making the value of a precomputed trajectory very 
questionable . 

Few studies have addressed themselves to the synthesis of a terminal-area 
trajectory. Although references 1 to 7 do not address this problem specifi- 
cally, they do provide insight into the terminal control process and are help- 
ful in identifying important problems. Reference 8 deals with the construc- 
tion of holding patterns and path-stretching maneuvers. Reference 9 is 
particularly relevant because it defines several of the problems studied here. 

This study investigates a technique for calculating the trajectory on 
board, with provision for recalculating it in real time whenever changes are 
needed. The practicality of this technique hinges primarily on the solution 
of two problems. The first is to establish what guidance problems occur most 
frequently in the terminal area and can be solved by an on-board system. The 
second is to develop efficient computational techniques for synthesizing the 
required trajectories. 
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The guidance problems to be solved are established from an analysis of 
air-traffic control procedures , instrioment flight maneuvers , and terminal-area 
route structures . They can be divided into three principal subproblems : a 

lateral (horizontal) maneuver^ a longitudinal (vertical) maneuver^ and a speed 
change. An important assumption made in this study is that these three prob- 
lems can be treated essentially independent of each other. 

A number of techniques are used to synthesize the trajectories; some are 
based on minimizing a specific performance function 5 others are derived from 
geometric considerations. In general, horizontal trajectories consist of con- 
nected segments of straight lines and circular arcs, while vertical trajec- 
tories consist only of connected segments of straight lines. Speed histories 
consist of sections of constant acceleration, constant deceleration, and con- 
stant speed flight. The choice of simple geometric shapes as basic building 
blocks yields trajectories that are easy to synthesize and that comply with 
all important aircraft maneuvering and performance constraints . There seems 
to be little reason to consider more complicated trajectories, at least for 
terminal-area operation of conventional aircraft. 

The basic information required to fly the synthesized trajectories is 
position, altitude, ground speed, and time. It is assumed that this informa- 
tion is available on board the aircraft. Although the description of the 
guidance laws derived herein is occasionally somewhat involved, programming 
experience has shown that the required computations can easily be performed on 
currently available airborne digital computers. (The development of a control 
system to fly the aircraft along the synthesized trajectory is closely related 
to standard autopilot design, and is not considered in this report.) 

DEFINITION OF GUIDANCE PROBLEMS FOR TERMINAL-AREA OPERATIONS 


Before trajectories can be synthesized, the set of guidance problems that 
occur in terminal-area operations must be defined and formulated mathemati- 
cally. Such a set can never be considered complete because there is no limit 
to the number of possible problems. Therefore, the objective is to choose the 
set as small as possible, consistent with the requirement that the most com- 
monly encountered guidance problems are included in it or can be reformulated 
in terms of problems in it . 

Terminal-area guidance can be viewed as the simultaneous solution of 
three types of problems: horizontal (lateral) guidance, airspeed management, 

and vertical guidance. To simplify analysis, this study assumes that each 
type of problem can be solved essentially independent of the other two. 
Although this assumption is not always Justified, it is motivated by observing 
that in commonly encountered flight conditions pilots tend to fly aircraft in 
a manner that allows them to solve these three types of problems indepen- 
dently. Although analysis of terminal-area guidance problems is given sepa- 
rately for each type of problem, in practice, the three types of problems do 
interact. The method of dealing with these interactions is given when the 
problem is encountered. 
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Horizontal Guidance 

The most complex guidance problems in terminal-area operations for con- 
ventional transport aircraft involve horizontal maneuvers. The origin and 
identity of these problems are most easily explained by a brief discussion of 
the route structure and air-traffic control procedures for an approach into a 
typical terminal area, that of San Francisco International Airport (SFO). 

The approach routes into SFO are indicated in figure 1 by the heavy 
broken lines. The most important feature is that they consist chiefly of con- 
nected sections of strai^t lines that terminate on the final approach. 
Transition from one straight-line section to another is achieved by means of 
smoothly rounded corners. The route structure can be specified by locating 
the corners defining each straight-line section and arrows indicating the 
direction of motion along each section. Corners often coincide with the loca- 
tion of a VOR station (such as Woodside) or at an intersection of radials from 
two adjacent VOR stations (such as Altamont). 

Under instrument flight conditions, air-traffic control closely monitors 
the movement of aircraft within the approach routes. The air-traffic control- 
ler sequences aircraft for landing, maintains separation between them within 
the routes, and vectors them onto the final approach with the desired spacing. 
In the existing ATC system, a controller normally observes aircraft on radar 
and issues the appropriate guidance instructions to pilots by radio telephone. 

Maintaining the proper spacing between aircraft is especially critical on 
the final approach and at points of confluence of two or more segments . Too 
small a spacing is unsafe and may disrupt the smooth flow of traffic, whereas 
too large a spacing results in a low -landing rate. Spacing distances for the 
final approach are typically in the range between three and seven miles, the 
choice depending somewhat on the individual controller as well as weather con- 
ditions . In weather conditions requiring instrument flight , controllers tend 
to space aircraft at greater distances than in clear weather. Spacing is 
achieved by controlling the speed of the aircraft or, if necessary, by vector- 
ing to control the length of the flight path to the point of confluence with 
the final approach. Vectoring the aircraft onto the final approach generally 
has the dual purpose of aiding the aircraft to acquire the ILS localizer and 
of ensuring the proper spacing between aircraft. Thus the route structure in 
the neighborhood of the point of confluence near the final approach (point A, 
figure l) is variable and not necessarily as shown; the shaded region approx- 
imates the area within which the route may vary. Having described the 
terminal-area route structure and the air-traffic control problem, the next 
step is to isolate the basic horizontal guidance problems that an automatic 
system would have to solve in guiding the aircraft during the approach. Five 
problems have been identified and are described in the following sections 
under the heading of problems (a), (b), (c), (d), and path stretching. 

Problem (a ) — Consider an aircraft that has just passed the corner at 
point B (fig. l) and is cleared for an ILS approach. (Assume that no other 
aircraft are on the approach routes.) The aircraft must be guided so that it 
will be on the final approach path and headed in the landing direction at 
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least by the time it passes the outer marker (OM). If it follows the approach 
route exactly as drawn in figure 1, it achieves heading alignment at point A, 
a few miles before the OM. In general ^ an aircraft must be in close heading 
alignment with the final approach at or before the OM, if the approach is to 
continue successfully. Thus, the general problem is to guide the aircraft 
from any initial position and heading to a specified final position and head- 
ing, with the final position and heading achieved -simultaneously. This prob- 
lem is illustrated by figure 2(a). 

Without some synthesis criterion there can be no unique solution to this 
problem, since for every combination of initial and final condition an infin- 
ite niimber of feasible trajectories can be found. Since a computer-oriented 
synthesis procedure requires imique solutions, the synthesis criterion was 
chosen as minimization of the path length of the trajectories connecting ini- 
tial and final points. This criterion was suggested by the shape of the 
approach route segments (fig. l) which show that once the trajectories have 
passed the last waypoint before the final approach they proceed by the most 
direct route to point A. 

An important constraint on the trajectories is that the bank angle 
required to fly along them not exceed some maximum value dictated primarily by 
passenger comfort considerations, typically 15° to 30°. The maximxim bank- 
angle constraint is, of course, equivalent to a constraint on the minimum 
turning radius; this equivalence is used in the synthesis procedure. 

Vvoblem (b ) — Major portions of terminal-area routes consist of straight- 
line sections; if an aircraft is already flying along such a section, a stan- 
dard lateral autopilot can keep it aligned. If, however, the aircraft is not 
on its assigned approach route, because the pilot or ATC have suddenly changed 
the route or a navigation error is discovered, then the problem arises of how 
to guide the aircraft onto the assigned route. This problem is the same as 
problem (a) when the new route must be entered at a specific point. Often, 
though, the only requirement is that the aircraft be flown onto the route as 
quickly as possible. Hence, problem (b) is defined as flying from any initial 
position and heading onto a straight line of specified location and heading. 

A solution to this problem is any trajectory that leads to the specified line 
without violating maneuvering constraints. As in problem (a), the selected 
trajectory should minimize path length without causing the aircraft to exceed 
a maximum bank angle. This problem is illustrated in figure 2(b), where the 
final line is assumed to be a VOR radial or the center line of the localizer. 

Problem (a ) — Another method of specifying an approach route is to choose 
a sequence of waypoints through which the aircraft is to fly consecutively. 

This method can also be used to generate holding patterns in a simple way. 

If trajectories connecting consecutive waypoints are chosen to minimize the 
path length, approach routes quite similar to those in figure 1 can be 
obtained by a Judicious placement of waypoints. 

Problem (c) is defined as flying from an arbitrary initial position and 
heading to a specified final position. As before, the synthesis criterion is 
chosen as the minimum path length, and the constraint is the maximum bank 
angle. This problem is similar to problem (a) except that the heading of the 
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aircraft when it arrives at the final position or waypoint is not specified. 
Figure 2(c) illustrates the problem for two waypoints. 

Problem (d ) — If the approach routes are specified by connected segments 
of straight lines as in figure 1, a procedure is needed for transitioning from 
one segment to another without excess overshoot or high bank angles. Problem 
(d) is defined as determining a control law for making a transition between 
two straight-line segments at any specified angle of intersection. An ideal 
transition between two lines is shown in figure 2(d). 

Path stvetohing — Air-traffic control was not considered in defining prob- 
lems (a) through (d), which arise in developing guidance laws for flying air- 
craft automatically along existing approach routes without regard to other 
aircraft. Air-traffic control requires knowledge of the location, heading, 
altitude, and speed of all aircraft operating in the terminal area. Since 
this knowledge is generally not available on board an aircraft, a comprehen- 
sive treatment of the air-traffic control problem from the point of view of an 
on-board system does not seem justified. It is reasonable, however, to con- 
sider on-board solution a portion of the air-traffic control problem that does 
not require explicit knowledge of other aircraft but would especially simplify 
the controller's job of spacing aircraft on final approach. This problem is 
to achieve precise arrival time of the aircraft at a designated point on the 
final approach or, in general, at any point on the approach route. Typically, 
such points might be chosen as the OM or at the junction of two or more 
approach routes (i.e., point A, fig. l). 

The preferred technique of achieving a specified time of arrival at a 
specified point on the approach route is through speed control, which does not 
require deviation from the assigned approach route. The use of speed control 
for this purpose will be discussed in the next section. However, speed con- 
trol alone is often inadequate, particularly if a change in the time of 
arrival is required close to the final point. The inadequacy of speed control 
is clearly shown when a speed lower than the stall speed is needed to achieve 
a delay in arrival time. Since late revisions in arrival time cannot be 
avoided, this limitation of speed control may be overcome by stretching the 
approach path as required to achieve the desired delay. 

Path stretching to delay arrival time is an essential element of all 
manual and automatic air-traffic control systems that have been studied in the 
literature (refs. 4 - 8 ). Although there is generally no unique trajectory 
that increases the path length by a specified amount, the selection of a par- 
ticular technique for path stretching is of great importance, since it is 
closely related to the method of air-traffic control. 

The effectiveness of path-stretching (and possibly also of path shorten- 
ing) techniques is often evaluated in terms of the timing accuracy they 
achieve in delivering an aircraft at the desired point. Of those evaluated 
in simulations, no clearly superior technique has so far emerged. For a 
description and comparison of five different techniques see reference 4 . 

Althou^ optimum methods of path stretching have not been developed, the 
method chosen should be able to increase the path length in a continuous 
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manner while preserving as much as possible of the approach route. Further- 
more, the calculations defining the stretched trajectory should be amenable to 
fast, on-line computer solution. Such considerations led to the particular 
method of path stretching described in detail later. 

To complete the definition of the problem, the path stretching is 
required to take place along a straight section of the approach on which a 
final point is specified. The path-stretching problem is then defined as 
selecting a trajectory that starts on the straight section, terminates on the 
final point, and has any specified arc length greater than the straight-line 
path. An example of a stretched path between two waypoints is shown in 
figure 2(e). 

Finally, the well-known holding pattern can be considered as an example 
of a stretched path. Generally, a path-stretching maneuver turns into a hold- 
ing pattern whenever the increase in path length specified is greater than the 
circumference of a holding pattern. Holding patterns are not discussed 
further in this report. 


Speed Control 

The primaiy purpose of speed control in the terminal area is to control 
time of arrival. In general terms, the problem is to calculate the speed-time 
or speed-position history of the aircraft so that it will arrive at a speci- 
fied point on the approach route at a specified time. This calculation 
assumes that the horizontal trajectory has already been determined by solving 
one or more of the trajectory problems defined in the preceding sections, so 
that the path length as well as the initial time and the final time tj^ 

are known. 

The speed-time histoiy is calculated subject to the aircraft’s opera- 
tional constraints on cruise speed V, acceleration, and deceleration. In 
addition to obeying these performance constraints, the initial speed Vq of 
the speed-time history must equal the current aircraft speed, and the final 
speed must equal a specified value. If the final point happens to be the 

OM, then may equal the final approach speed with fully lowered flaps. 

The "simplest" speed-time history that can solve this problem consists of 
three segments: an acceleration or deceleration segment to the cruise speed, 

a cruise-speed segment, and an acceleration or deceleration segment to the 
final speed, as in figure 2(f). 

Calculation of the speed-time history may reveal the need to apply path 
stretching to the horizontal trajectory. For example, the path length will 
need to be stretched if the calculated cruise speed is less than the minimiim 
speed. 


Vertical Guidance 

For the arrival operations considered in this report, vertical guidance 
is restricted to descent from cruise. (For departure operations, similar 
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considerations are applicable to ascent to cruise.) Here it is assumed that 
the horizontal trajectory and the speed-time history have been obtained as 
solutions to previously defined problems . 

A typical altitude history, illustrated in figure 2(g), consists of seg- 
ments of constant altitude and of descending flight. It commences at cruise 
altitude h^ and terminates at a final altitude hf approximately equal to 
the height above the runvay of the ILS glideslope beam center at the OM. 

The descent profile provides for maintaining the cruise altitude as long 
as possible, then descending to the final altitude at a specified descent rate 
and at a constant forward speed. If any speed change is needed near the final 
point, it should be made in level flight; the descent profile must be planned 
to allow for the required period of level flight after the final altitude is 
achieved. Thus, the speed-time history (and horizontal trajectory) must be 
calculated before the altitude profile. 


Summary 

A complete terminal-area trajectory, consisting of the position, speed, 
and altitude of the aircraft as a function of time, is obtained by solving 
sequentially the horizontal, airspeed, and vertical guidance problems that are 
defined in this section and summarized in table 1. Although these are not the 
only problems that can be defined, they are sufficient for calculating the 
most frequently occurring trajectories in present-day terminal-area opera- 
tions. (important additional problems could be defined if information on the 
locations and trajectories of other aircraft in the area were available on 
board the aircraft, but such problems are not discussed in this study.) 

Solutions to the various problems discussed in this section are developed 
without explicitly considering the wind. Thus, they yield values for ground 
heading, ground speed, and ground acceleration time histories. If the wind 
speed and direction are known, as we assume in this study, these quantities 
can easily be converted to aircraft heading, airspeed, and acceleration with 
respect to the air mass. 


TRA.JECTORY SYNTHESIS FOR PROBLEM (a) 


In reference 9, trajectories were derived for problem (a) which are of 
minimum path length and are constrained by the minimum turning radius . \>Jhen 
the initial and final points are separated by less than four turning radii, 
such trajectories can contain up to three partial turns of minimum turning 
radius. When the initial and final points are separated by more than four 
turning radii, the trajectories need only contain a straight-line segment with 
a partial turn of minimum turning radius at each end. In this section, a 
computer-oriented algorithm for calculating the trajectory is developed based 
on certain simplifying assumptions and on a geometric interpretation of the 
guidance problem. 
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The trajectory patterns to be used in the development of the algorithm 
consist of at most tvo circular arcs and a straight-line segment* The three- 
circular-arc trajectories of reference 9 niay occur in the minimum-arc-length 
solution of problem (a) if the distance of separation between initial and 
final points is less then four turning radii. However 5 they are excluded here 
since they would increase the complexity of the algorithm and are expected to 
occur infrequently in practice. 

Another difference between reference 9 and the present development con- 
cerns the optimality of the trajectories generated by the algorithm. Refer- 
ence 9 was devoted to a mathematical derivation of the minimum arc problem for 
problem (a). Here, for computational simplicity, the trajectories are not 
selected to minimize, in the strict mathematical sense, a specific performance 
function such as the arc length. Rather, the key criterion for selecting tra- 
jectories is geometric constructability; that is, the feasibility of con- 
structing a particular type of trajectory for given initial and final condi- 
tions. Because the trajectories for which a feasible construction is sought 
are known to be optimum under some conditions, as shown in reference 9 5 the 
procedure can be expected to yield an efficient trajectory, though not neces- 
sarily the true optimum one in. the minimum arc sense. The disadvantage of not 
always generating trajectories that minimize a performance function is out- 
weighed by the geometric insight which the algorithm casts upon the horizontal 
guidance problem. 

It was stated previously that the trajectories should obey a constraint 
on the maximum bank angle. At constant airspeed V a constraint on the mini- 
mum turning radius is equivalent to a constraint on the maximum bank 

angle 4 tiax* relationship between these constraints is given by the 

following equation: 


R . = 1) 

min g tan 6 

^ ^max 

where g is the acceleration of gravity. The value of V used in equation 
( 1 ) should be the maximum ground speed encountered during the maneuver. An 
acceptable estimate would be to choose it as the sum of the maximum airspeed 
in the terminal area and the wind speed. In further calculations, the turning 
radius is treated as a parameter whose value is arbitrary and is related to 
the bank-angle constraint. 

The trajectories in this study are chosen from the 13 basic patterns 
shown in the left column of table 2. (The additional patterns shown in the 
right column of table 2 are required for the complete minimum arc solution of 
problem (a), and are not used in the algorithm.) The patterns used are com- 
posed of three or fewer segments ; each segment consists of a portion of either 
a straight line or a circle. Under the maneuver symbol column in table 2, 
each pattern is assigned a maneuver symbol, characterizing the consecutive 
maneuvers required to fly the trajectory. For instance, L S L stands for 
left turn, straight flight, and left turn. These symbols are used throughout 
the text to refer to a particular pattern. 
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Pattern Parameters 


The assignment of a flight pattern depends on the initial and final con- 
ditions of the problem, represented hy the coordinates of the initial point 
and the initial heading and by the coordinates of the final point and the 
final heading. For the purpose of pattern assignment and the calculation of 
pattern constants, the initial and final conditions can be reduced to three 
parameters or states. 

The symbols Pq and P^. denote the initial and the final point and 
and ijjf the initial and the final heading, respectively. Since is 

normally the runway heading, the initial heading is conveniently measured with 
reference to ipf. The initial heading thus measured is denoted by i(; and is 
given by - i|;f. Unit vectors n^ and n^ have directions and 

and emanate from Pq and P^, respectively (fig. 3). If n and rif are not 
parallel or antiparallel, lines drawn through Uq and fi^, termed the "initial” 
and "final” lines, will intersect at a point P^. The case where these lines 
do not intersect is referred to as the singular case and is considered 
separately elsewhere. 

The pattern parameters for the nonsingular case are ip , d , and d^. A 
positive ip is defined by rotating n^ clockwise into The range of \p 

is -7T < ^ 7T, with ip in radians . The other two parameters are dQ and dr 

where dQ is the signed distance Pq^v is the signed distance P^P|». 

The sign of a distance is positive if the corresponding unit vector points 
away from P^. The geometric construction of these parameters and a typical 
trajectory are illustrated in figure 3. 

The two parameters dQ and d^ are not directly measured, but must be 
computed from available navigation measurements. The computation is aided by 
relating them to a commonly used Cartesian coordinate system. As illustrated 
in figure 3, the origin of the coordinate system is at P|»; its x-axis points 
in the direction of nf, and its y-axis points to the right of the reference 
direction. This coordinate system is often used in aircraft guidance and con- 
trol studies. Let (x,y) be the coordinate vector of Pq in this system, nor- 
malized by using the turning radius R as the unit of distance: 

(x*,y*) = (x/R,y/R). Then dQ = y*/sin , and d^ = -x’ + y*/tan ip, or 


d 


o 


y . 

R sin 


( 2 ) 


d 


f 



y . ■._ 

R tan 


(3) 


The quantities and d^ are also normalized distances. Equations (2) and 
(3) are valid for all ^ except = 0 or = tt , corresponding to the singu- 
lar cases. In an actual computer mechanization, ip will be considered zero or 
TT when it lies in the neighborhood of these values. 
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Switching Points 


Our objective in the next several sections is to describe how flight pat- 
terns are assigned as a function of the three parameters d^, and d^ and 
to determine those combinations of parameters, called the switching points, at 
which changes from one type of pattern to another occur. Although the switch- 
ing points can be functions of all three parameters, the pattern assignment 
problem is simplified by choosing many of the switching points as dependent 
only on ij;, and the others as simple functions of and d^. The pattern 
assignment problems are different for nonsingular and for singular values of 
ijj; the assignment problem for the nonsingular case is considered first. 

All switching points are defined by points of tangency involving the ini- 
tial or the final line and certain circles which arise in the geometric con- 
struction of flight patterns. The six switching points on the initial line 

are P-.P ^.P-,P^,P+, and P +, and the six switching points on the 
oa o 6 oa og oa oa 

final line are ^fa’ ^f6’ ^fJ’ ^fa’ ^fy’ directed distances measured 

from the origin at P^ to these switching points are 9 ^06’ ^oa’ ^06’ 

d +, d + for those lying on the initial line and d - , d -, d d^+, d_+, 
oa oa ^ 

d^+ for those lying on the final line. The "initial" circle is the circle 

tangent to the initial line at Pq , and the "final" circle is the circle 
tangent to the final line at P^. 


Flight Patterns for Normal Maneuvers 

The "normal" flight patterns apply, and the derivation of switching 
points is simplest, if the initial and final points are far apart compared to 
the turning radius. Figure k illustrates the trajectories for this case with 
an acute and an obtuse value of ^ for 0 ^ ip it. Each part of the figure 
shows four trajectories corresponding to four pairs of initial and final con- 
ditions, with designations corresponding to those in table 2. (Circles 
involved in the construction of the trajectories are completed by broken 
lines.) The region of large separation between initial and final points can 
now be defined more precisely as the set of initial and final conditions for 
which neither of the two initial circles intersects either of the two final 
circles. In this region, the four trajectories shown in figures h(a) and U(b) 
can always be drawn. The boundaries of these regions can be defined by study- 
ing the figures as either the initial or final point moves along its line. 

The trajectoiy starting at Pq^ and ending at P^ ^ is of the LSL 
type. As Pf 1 moves down, the trajectory changes to the RSL type when the 
left initial circle becomes tangent to at switching point ^h^ dis- 

tance of P + from P is d_+ = tan ip/ 2 , 
fa V fa 

The same trajectory, as P^^ moves to the right, changes from LSL to 

LSR when the left initial circle is tangent to the line d^ at P P - 

° i oa oa 
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and are both on the same circle. The distance of P - from P-- is 

fa oa ^ 

d - = -tan i|;/2. 

The final two switching points are defined from geometric convenience in 
contrast to the first two, where the switching must take place at the defined 
points. A trajectory starting at P^^ and ending at P^^ is of type RSL. As 

P^^ moves to the right, the right initial circle first becomes tangent to the 

final lines at P _ , at d ^ = -l/tan(ijj/2) . 

Op Op 

The last switching point, modifies the trajectory connecting P^^ 

and moves toward ^^ 2 » switch from type LSR to type RSR 

occurs at with this switching point is not a geometrical 

necessity, since both LSR and RSR trajectories are usable for all points on 
df. However, the switch from LSR to RSR at some point near P^ yields tra- 
jectories of shorter path length. Since P^+ is a geometrically significant 

point for near-distance maneuvers, it was chosen for this purpose also for 

minimizing the total ntimber of switching points. Its location on dj^ is 

established by locating P^^ so that the left initial circle of P ^2 a-nd the 

right final circle of P^^ are tangent at only one point as P^^ moves 

toward P^^ . P „ is then at P + , and P^, is at P.+ (fig. 4}, To con- 
fa 02 oa’ fl fy ^ 

struct the circles establishing these points, draw a circle with center at the 

intersection of two lines, one line three units to the right of d^ and the 

other line one unit above d^, tangent to d^ at second circle 

can now be drawn tangent to this circle, and also tangent to d^ at 

Thus, d^+ = 3/tan ij; + l/sin ij; and d^+ = 3/sin 4^ + 1/tan ij;. 

From symmetry, trajectories for the negative range, -tt < i); < 0 , can 
be taken as the mirror images of trajectories for the positive range 
reflected about the d^ axis. The reflection of trajectories for negative 
values is essentially equivalent to the introduction of another switching 
point at = 0. This reflection procedure for negative ifj is also valid for 
near-distance maneuvers . 


Fli^t Patterns for Near -Distance Maneuvers 

If the initial and the final points are so close that at least one pair 
of initial and final circles intersect, the flight patterns given in the pre- 
ceding section may become invalid. Flight patterns for this condition are 
referred to as near-distance maneuvers. Near-distance maneuvers rarely occur 
in terminal-area operation of aircraft and should always be avoided since the 
resulting disruption of aircraft on landing approach can cause landing delays. 
Nevertheless, events that generally require near-distance maneuvers, such as 
go-arounds and emergencies , do occur in a small but not negligible percentage 
of takeoffs and approaches. Therefore, no automatic terminal guidance system 
is conqplete unless it can also handle such events. 
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Pattern selection rules for near-distance maneuvers are somewhat 
involved. Readers not interested in following these details may go directly 
to the next section, on switching diagrams , where the selection rules are 
summarized. 


When Pq is to the left of ‘t'wo groups of switching points are 

needed to define the trajectories. The first group contains the points 
P^|, and P^-. If in figure k{a) or U(h) the points P^^ and P^^ move simul- 
taneously toward the origin, they will reach positions where the straight-line 
segment of the RSL trajectory vanishes and the right initial circle becomes 
tangent to the left final circle. There are generally two possible final 
circles tangent to the line d^ and to the initial circle; in figure 5(a), 
points of tangency of the final circles with d^ are denoted by P^- for the 

further of the two points and by P^^ for the nearer of the two from the 

origin. The positions of these two points are functions of the position of 

P ; as P moves back to the left to P -, the points of tangency P«- and 
ol ol oa fy 

P^^ merge into a single point. For initial points to the left of ^^”5 ini- 
tial and final circles do not intersect, and the region is as discussed in the 
preceding section. The upper limit of is P^+ , reached when P^^ is at 

P -. 
oa 


The second group contains the points P^- and P 


■ o6 ' 


Point P - is 
fa 


defined by a circle tangent to the negative d^ axis and the positive d^ 
axis, as illustrated in figure 5(a). The definition of P^^ depends on the 

value of !(;. For 0 < < tt/2, P^^ is obtained by moving P^^ toward P^- 

\intil P^J coincides with figure 5(c). For tt/2 < ij; < tt, P^^ is 

obtained similarly by moving P^^ toward P^- until P^- (instead of P^;j) 

coincides with P^-, in figure 5(d). This definition of P . is continuous 
fa oo 

at = 7r/2, the angle common to both definitions. Thus, for small P 


o5 


lies near P -. As ip increases 
oa 


■ o6 


moves toward P coinciding with 
oa 


p _ at ip = 7t/2. For lii > 7t/2, P . reverses its direction and moves in the 
oa r / ^ ^ oo 

positive direction. 


For final points in the trajectory assignment depends on 

values of and on the initial point. If 0 < < 7r/2, and the initial point 

is in [P -,P .], an RSR trajectory is used as shown in figure 5(a). For ini- 

tial points in trajectory is used for final points below 

P - as in figure 5(a), and an LSR trajectory is used for final points above 
fa 


Pf-, as in figure 6(a). If tt/2 < < tt, an LSR trajectory is assigned for 

Pq in » either an LSR or an RSR type for P^ in [P^^,P^-], an LSR 



trajectory is used for final points above on RSR type for final 

points below 

The effect of on pattern assignment is to modify the definition of 

P^-. Whenever Pq is at or to the right of corresponding P^- is 

at negative infinity on the d^:. axis. If P* is to the left of P q, P«- 

^ o p 1 Y 

is located by the procedure illustrated in figure 5(a). 

For initial points in final points not in tra- 

jectories are assigned so as to be compatible with those in the preceding sec- 
tion. Thus j if P^ is in type is used; if P^ is above 

P^+, an LSL type is used. Final points can be below P^- only if P^ is in 

[P -,P since otherwise P_- is at negative infinity. In addition, P „ 
oa’ o6 fy o3 

and P^- interchange the order of their positions on the d^ axis as 

passes through tt/ 2 radians. However, for all values of i(^, the RSL type is 

used if Pp, is in [P -,P and P« is below P«-. This completes the tra- 
° oa’ oB f fy 

jectory assignment for the set of initial points to the left of P - and all 

points on d^. 


Trajectories for points to the right of 


P - involve switching points 
oa 


P^+, P +, P +, and arise as a point on 
fy’ oa’ oa’ ^ 

moved toward 


d^ , far to the right of the origin, is 


P In figure 4, as P ^ and P^, move simultaneously toward 
oa ^ ’ 02 fl 


the origin, they reach positions where the initial and the final circles 
become tangent. As discussed previously, past this point there are two final 
circles tangent to the initial circle. However, unlike the case discussed 
earlier, only the tangent point farthest from the origin is of interest. This 
point (P^+) "tbe trajectory that defines it are illustrated in figure 5(b). 

The position of P^^ , which yields only a single tangent point between initial 

and final circles, as P^^ and its right final circle are moved along d^, is 

defined (fig. U) as P +. 

oa 


For P to the left of P + , P^+ is a function of 
02 oa fy 

both ^ and whereas to the right of PqJ> P^+ is arbitrarily assigned 

the value corresponding to P^^ at P^-* For all initial points to the right 
of P^- and final points above P^:^» s-n LSR type is chosen, which is simply 
the continuation of the type used in figure k between P_ and P^, . Simi- 


02 


fl ' 


larly, an RSR type is used for all initial points to the right of P^+ and 
final points below This type is shown in figure 5(b) as the trajectory 


connecting points P^^ and P^^ “ switching point PqJj which is the com- 

panion point to P - and is illustrated in figure 5(b), will be discussed in 
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the next paragraph. ) The RSR type is also used for final points below P - 

if the initial points are in [P -,P +]. 

^ oa’ oa 

The trajectory assignment problem is concluded by specifying the trajec- 
tories for Pq in an<i Pf in ^ switching 

point at P^+ is shown when P^^ in figure 5(b) is moved into coincidence 

with P+ while maintaining P„„ between P_- and P_+. At the point of coin- 
oo 12 icy ly 

cidence, the RSR trajectory assigned to points to the right of P + degene- 
rates into a RS trajectory, thus establishing P^+ as a limiting point of the 
RSR type. Similarly, the companion point to Pq+> is a limiting point of 

the RSR type. If, in figure 5(b), P^^ is in and P^^ is below P^-, 

the RSR trajectory assigned in the preceding paragraph to this combination of 
initial and final points degenerates into an SR type as P^^ moves into coin- 
cidence with Pf 2 sli^tly above P^-» and the left initial 

circle and the right final circle do not intersect, an LSR type is used, shown 
in figure 6(b) between P^ and P^. The final switching point, denoted by 
P^^, then occurs as P^ is moved up \mtil the left initial circle and the 

right final circle become tangent. For final points above "the LSL type 

in figure 6(c) is used. The range of P^^ is "the lower limit 

being reached as P^ approaches P + and the upper limit as P approaches 

P -. 
oa 


As before, these rules must be slightly qualified to taJce into account 

the effect of P ^5 which moves to the right of P - for ib > 7t/2. If P 
og^ oa o 

lies to the right of P , no changes in the rules are needed. For P^ in 

Op o 

[Po5»Poe^’ P^~ is located with respect to P^. If P^ is above 

is disregarded and the established assignment rules are used. If is 

below the previously established RSL trajectory is used. This concludes 

the assignment of trajectories for all combinations of the pattern parameters. 


Derivation of Switching Functions 

The next step is to derive the analytical expressions not yet given for 
some switching points. An essential step in the derivation is to solve for 
various parts of the right triangle P^, Cq, P^, shown in figure J, Since all 
circles in the construction have unit radius , = /l + d^ and 

11^2 = ^ + tan“^ 1/^0 where d^ is negative in this example. Now 


CP = y/1 + d^ sin ^2 j 
o n o ’ 


l6 



P P = /l + dZ cos \p2 
n V o ^ 


' VfJ = V* ■ 


These expressions suffice to calculate and 

d^J = d^ cos ip - sin ip + A - (d^ sin ip + cos ip + l)^ 


d^- = d cos xp - sin ip - A - ( d sin \p + cos ip + l)^ 
f Y o ^ r ' o 


(h) 

( 5 ) 


The expression for d^- can he obtained by equating to zero the expres- 
sion for solving the resulting equation for d^, and substituting d^- 

for d : 
o 


oa sin Ip tan ip 


( 6 ) 


The expressions for d + and d - are: 
^ oa fa 


d + = tan — 
oa 2 


d - = -tan 
fa 2 


( 7 ) 

(8) 


If 0 < < Tr/2, d^^ is obtained by equating expressions for d^^ and 

d^- and solving the resulting equation for d^, which is then interpreted as 


o6* 


^o6 ~ 2 - sin 


-v-( 


1 + cos Ip - tan ^ sin ij;] ^ 


( 9 ) 


Similarly, for tt/ 2 < t/; < tt , d^^ is obtained by equating d^- and How- 

ever, after the proper sign is chosen in front of the square root expression 
occurring in the solution of the quadratic equation, the equation obtained is 
the same as equation (9)* 

The expressions for d^+ and d^^ can be derived in the same way as d^- 
and They may also be deduced directly from d^J and d^- by symmetry: 


d^+ = d^ cos ip + sin ip + A - (-d^ sin Ip + cos i|j + 1)2 (10) 


IT 



(11) 


= d cos 
f 6 o 


}p + sin \p ~ A"- (-d sin + cos \jj 

n 


A complete collection of all switching functions, along with their 
domains of definition, is given in table 3- 


Switching Diagrams 

The various regions of trajectory patterns in the three-dimensional 
parameter space (i|), dQ, dp) can be illustrated with switching diagrams which 
consist of planar cuts through this space obtained by holding one of the three 
parameters constant. When ip is constant, the switching diagrams are as 
shown in figures 8(a), (b), and (c) for three values of ip: tt/U, tt/2, and 

3tt/ 4. Each region is labeled with the pattern type established for it in the 
preceding sections. On the region boundaries, the four types of three- 
segment patterns often degenerate to one- or two-segment patterns, as indi- 
cated by the coded lines. In addition to giving a compact representation of 
pattern regions, the switching diagrams find use in developing decision logic 
for pattern type selection. 

The switching diagrams show that for \p = n/h and 3 tt/U, d^^ = The 

line segment d^- is terminated on the left by d^- for < 3'ir/U and by d^^ 

for > 3 tt/ 4, as shown by the two inserts in figure 8(c). No other changes 

in the structure of the switching diagrams occur as is varied in the 

range 0 < ij) < tt. 

Flight Patterns for if; = 0 and i() = tt Radians 

For the two singular values of ip^ \p = 0 and if) = tt radians, the param- 
eters do and dp are not defined. Instead, the parameters used are the rec- 
tangular coordinates of the initial point, P^, expressed in the X,Y coordi- 
nate system of figure 3. This coordinate system is centered at Pp and has 
its x-axis pointing in the direction of the final heading. Unlike the param- 
eter space (do, (ip), the parameter space (X, Y) preserves the geometric shapes 
of figures under mappings; therefore pattern assignment and switching diagrams 
can be combined in one figure. 

Switching diagrams, as well as example trajectories, are given for 
if) = 0 in figure 9(a) and for if) = tt radians in figure 9(b). For both values 
of if), the switching curves which separate the parameter space into regions of 
different patterns are either straight lines or portions of circles with 
radius 2. For if) = 0, one circle is centered at (0,-2), the other at (0,2). 
For if) = 77, the one circle is centered at the origin. These circles define 
the locus of initial points for which an initial circle is tangent to a final 
circle. The switching diagrams are symmetric about the x-axis. This sym- 
metry is also present in the pattern assignment; for a given value of X, tra- 
jectories in the negative Y plane are the mirror images of trajectories in 
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the positive Y plane* Several sample tra,1ectories are superimposed on the 
switching diagrams to illustrate the pattern assignment. 


Derivation of Path-Length Equation 


In preceding sections , rules were derived for choosing a particular tra- 
jectory for each initial position and heading of an aircraft. The next step 
is to derive equations for calculating the path length and the heading changes 
in the first and final turns for each trajectory. For the RSR trajectory 
between and in figure 10(a), P^ is at (x,y), the initial heading is 
and the turning radius is R- Then: 


P C = P^P + P C 
f or f o o or 


X + R cos 


+ |-j, Y + R 


( 12 ) 


Further, the vector C is the difference between P«C and P_C : 

^ or f r T T T' T 


f fr 


f or 


C c, 

or fr 


-X - R cos(t|) + |-j, R - R sinf^' 


(13) 


The heading of the vector denoted by , equals the heading of the 

straight-line segment, 4)g , whereas the length of this vector equals the length 
of the straight-line segment : 




s 


heading of straight- 
line segment 


, -1 R - R sinCit) + (tt/2)] - Y 

-X - R cos[t|) + (tt/2)] 


(lU) 


P1P2 


length of straight- 
line segment 


/{X-^R cos [(//+( tt/2) ] }^+{R-R sin [ j|/+( tt/2 ) ]-Y}^ 

(15) 


The arctangent is evaluated in the conventional manner as an angle in the 

range 0 < ^ < 27 t, The heading change in the first turn is seen to be 


ij;. = 4) - 

1 S 

4 ^ if 4»^ > 4< 


( 16 ) 

4). = 2tt + 
1 

'f’s - 

if 4' < 

s 


(IT) 

The path length in the first turn is 
4^^ and the path length of the final 

then PqPi = 4^il^ • 
turn are 

Finally , 

the angle 

4/^ = 2tt - if 

4^ >0 

s 

P2Pf = 

R4^f 

( 18 ) 

•H 

0 

II 

4j = 0 
s 

PaPf = 

0 

( 19 ) 
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The trajectjDiy constants for an RSL trajectory, as in figure 10(h), begin 

( 20 ) 


with vector C : 

or fe 


“^or^fe ^f'^or 


which can he written in Cartesian coordinates: 


C C, 
or fe 


-X - R cos(i|; + Ij, -R - R sinU + |- ) - Y 


( 21 ) 


For the trajectory to exist we must have H ~ existence of the 

trajectory is guaranteed if the pattern selection rules established in the 

preceding section are applied. *^ e next ste p is, t o calcula te ^ ^ the angle 

between the vectors PTro and C ; since C P = (l/2)||C || , 

^ ^ or fe ’ or c "or fe ’ 


\p^ = arc sin 


I ' Ofe 


( 22 ) 


The length of the straight-line segment and its heading are therefore 


P 1 P 2 


length of straight- _ 2R 
line segment tan 


(23) 


, _ heading of straight- , . , , 

^s line segment ^ 


if + ij)’ < 2 t7 

= - 2tt if 1(^1 + > 2ir 


(2h) 


where is the heading of the vector 

-R - R sin[i|j + ( 77 / 2 )] - Y 
♦1 ' -X - E cost* + U/2)J 


( 25 ) 


Finally, the heading changes and arc length of the initial and the final turns 
are : 


^i = 


'Ps - 'p 


+ 2-rr - Ip 


\p > ip 
s 


, % 'P 


PqPi = 71;. R 


( 26 ) 


’^f = '^s ’ 


P2Pf - \p ^ 


( 27 ) 
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For the LSL and LSR patterns, the same techniques used for the RSR pattern 
above are applicable* Table 4 summarizes path-length and heading change equa- 
tions for the four major pattern types. When possible, the equations are 
simplified by application of trigonometric identities. Patterns with fewer 
than three segments are special cases of the three-segment patterns and are 
not included in the table. 

The major computational units in the horizontal guidance algorithm are 
summarized in the flow chart given in figure 11. Reference in the chart to a 
figure or table in the text identifies a major unit and implies that logic 
must be developed to perform the indicated function. For example, entering a 
switching diagram to select a pattern requires the development of a flow chart 
of inequality tests performed on the switching functions given in table 3 to 
locate the correct pattern region. Examples of such flow charts for the 
switching diagrams of problems (b) and (c) are in the next section. 


TRAJECTORY SYNTHESIS FOR PROBLEMS (b), (c), (d) 


Optimum control laws for the minimum-path-length trajectories of prob- 
lems (b) and (c) are derived in reference 9. To complete the solution of 
these horizontal trajectory problems, equations are required for the path 
length and the turning angles of the optimum trajectories. In addition, com- 
puter algorithms are needed for the implementation of the control laws. Fur- 
ther, a solution to problem (d) is necessary. As in the previous section, 
distances are normalized to units of turning radius. 


Trajectory Synthesis for Problem (b) 

Problem (b) was to fly from any initial position and heading to intercept 
and fly along a line of specified heading. The mi nimim- arc -length control law 
for this problem, from reference 9, is shown as the switching diagram in fig- 
ure 12. Coordinates are d and ; d is the perpendicular distance in units 
of turning radii from the final line, and is the heading measured with 

respect to the heading of the final line. The coordinate d is positive for 
initial points to the right of the final line if one faces in the direction of 
the final heading. The heading ip is positive for a clockwise rotation from 
the direction of the final line, and its range is -tt < i/; < tt. Solid switching 
lines with arrows indicate the direction a phase point travels toward the ori- 
gin. The switching lines and the dashed lines separate the phase plane into 
various regions. Each switching line and region is labeled with its respec- 
tive pattern type. Reference 9 shows that the optimum trajectories for this 
problem consist at most of three segments , with the middle segment a straight 
line and the first and last segments portions of minimum-radius circles . 

In the phase plane, right and left turns with imity turning radius are 
parameterized by the equations 

^o ^ ^ (right turn) ( 28 ) 
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(29) 


lllllllllllli 


d. (if/) = d - cos \p + cos if/ (left tiorn) 

36 O O 

where do and are the phase points where the tuirns are initiated. 
Straight-line flight generates a line parallel to the d axis. The switching 
lines are special cases of these phase plane trajectories and are labeled Sj 
through Ss in the switching diagram. The equations for the switching lines, 
along with their domain of definition, are required hy the computer algorithm 
to select the correct pattern type for a given initial condition: 


Sj (i|/) 

= 1 - COS Ip 

Ip ^ [ 0 , it/ 2 ] 

( 30 ) 

S 2 ( 4 >) 

= 1 + cos \p 

-tt/2 £ rp ^ tt 

(31) 

Ssi'P) 

= -1 + cos Ip 

\p 1 [0, -tt/2] 

(32) 

34(d) 

= -i\l 2 

d > 1 

(33) 

85(d) 

= 7t/2 

d < -1 

(34) 

Se(')') 

= -1 - cos Ip 

-TT < \p < tt/2 

(35) 


The dashed lines which also form portions of the boundaries between pattern 
regions do not represent phase-plane trajectories. Instead, they define the 
set of initial points for which two trajectory types yield equal path length. 
For example, the RSR trajectory and the LSR trajectory yield the same value 
for the path length for initial conditions on the dashed line segment 
4) = tt/2, d > 1. For definiteness, one of the minimum-path-length trajectories 
is arbitrarily assigned to such points of ambiguity. The phase trajectories 
for eight initial conditions (P^ - Pj^) are traced out in figure 12, and their 
corresponding geometric patterns are shown in figure 13. 

Parametric equations (28) and (29) for the phase-plane trajectories, 
together with the switching diagram, are used to derive equations for turning 
angles and for the length of the straight-line portion of the various pat- 
terns . Table 5 lists these equations and their region of validity for all 
trajectory types beginning with a left turn. Equations for the mirror-image 
trajectories are obtained by reversing the signs of ip and d. For example, 
if d and ip are such that the switching diagram specifies an RSR trajectory, 
the signs of d and \p are reversed, then entered into the equations for the 
LSL trajectory given in table 5- 

Finally, the switching diagram (fig. 12), equations (30) through (35), 
and the equations in table 5 can be used to construct a flow chart for trajec- 
tory synthesis (fig. 1^). The input to the flow chart is d and ip ^ and its 
output is the trajectory type and parameters. The first part of the flow 
chart consists of a sequence of inequality tests on d and \p to determine the 
pattern type, and the second part consists of the calculation of the trajec- 
tory parameters. The flow chart logic was simplified by making use of the odd 
symmetry of the switching diagram. If the given d and ^ require a pattern 
beginning with a right turn, the algorithm reverses their signs, identifying 
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the situation by setting index = 2, After the pattern ■oarameters are com- 
puted, i;^ is checked; if it is 2, the pattern is reversed to its mirror image 

Trajectory Synthesis for Problem (c) 

Problem (c) is to fly from a given initial position and heading to a 
specified final point with an imspecified final heading. The control law that 
minimizes arc length, derived in reference 9, is shown in figure 15(a), with 
distances normalized to the unit turning radius. The location of the final 
point is expressed in a coordinate system which has its origin at the initial 
position, Pq 5 its x-axis pointing in the direction o^^ the initial heading, 
and its y-axis pointing to the right. The switching lines in this control 
law are the x-axis and two circles with unit radii and centers at (0,-l) and 
(0,l). Trajectories for points on the negative x-axis can be taken as either 
LS or RS with no change in path length; they are arbitrarily made RS. Example 
trajectories for four final points are shown in figure 15(b). 

Equations for the path length and turning angles , and the flow chart for 
trajectory synthesis, are determined next. Symmetry in the control law 
reduces the complexity of the trajectory constants and algorithm. Only the 
positive y-plane need be considered, since trajectories for the negative 
y-plane can be obtained by reflection about the x-axis. Then the switching 
logic need only determine whether a final point is inside or outside the unit 
circle at (0,l), and the path length and turning angle equations are easily 
derived by simple trigonometry. Table 6 summarizes the equations and their 
region of validity; denotes the angle in the initial turn, the angle 

in the final turn, and the length of the straight-line segment. (The arc- 

tangent is an angle \|j in the range -tt < ^ < tt.) The flow chart for trajec- 
tory synthesis, shown in figure l6, is an implementation of figure 15 (a) and 
table 6. It has as inputs the x,y coordinates of the final point, and as 
outputs the trajectory pattern, turning angles, and path length. ^or purposes 
of developing this flow chart, each pattern was assigned an integer, which is 
represented in the chart by the index i-p. The correspondence between the 
pattern types and this index is given in the chart. 


Trajectory Synthesis for Problem (d) 


Problem (d) is to make a smooth transition from flight along one straight 
line segment to flight along a second, intersecting straight-line segment. 

The distance from the aircraft to the point of intersection along the first 
segment is Xj^ , and the distance from the point of intersection to a desired 
final point on the second segment is x^^. Angle is the heading difference 

between the two segments; since the segments intersect, ifj is neither zero nor 
7T. To make the transition without overshoot, the aircraft should begin a cir- 
cular turn of radius R at a distance before the intersection point given by 


x^ = R tan 


M 

2 


< ij; < 7T , or 
-7T < < 0 


( 36 ) 
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When proper account is taken of the circular segment, the total path length 
between and Xf is 




T 


= X. 
1 


+ 


- 2R tan 


M+M 

2 R 


(37) 


CONTROLLED TIME OF ARRIVAL AND VERTICAL GUIDANCE 


The problem of achieving precise arrival time of an aircraft at some 
point on the final approach or any other specified point of the approach route 
generally consists of selecting both a speed profile and a horizontal trajec- 
tory of appropriate length such that the aircraft', starting with initial speed 
Vq arrives at the end point with speed Yf in the time interval t. As 
explained earlier, we first attempt to achieve the desired arrival time with 
speed control alone, assuming that the horizontal trajectory, and therefore 
the path length, is fixed. If it is fo\md that speed control cannot achieve 
the desired arrival time because the aircraft's performance capabilities are 
exceeded, the path length will be modified as required to control the arrival 
time. Solution of the controlled time of arrival problem therefore involves 
two algorithms, one to select the speed profile, ass-uming a fixed length of 
trajectory to the final point, the other to change the length of the trajec- 
tory by an appropriate amount. 


Speed Control 


The problem of speed profile selection is considered first. The time to 
the final point and the length of the path are denoted by t and L, respec- 
tively. In addition to beginning with Vq and ending with Vp , the speed V 
along the profile must at all times obey the constraint Vmin ^ Y < V^ax* 

Speed changes are made by constant acceleration Aa (nximerically positive) or 
constant deceleration A<j (mmerically negative). These six quantities, t, L, 
V-tni n < Vmavi , and A(j, are treated as parameters in developing logic for 
speed profile selection. The speed profile consists of a maximum of three 
segments: an acceleration or deceleration segmeftt starting at Vq, a constant 

speed segment, and another acceleration or deceleration segment ending at Vp. 
Unless one or more of the segments has zero time duration, the three segments 
can be arranged to yield four possible speed profiles. The profile to be 
used is selected by comparing the path length L with two distances, L^ and 
L2 : 


Li 



+ 


2A= 

2^d 


""o ^""f 




(38) 

(39) 


2k 



r 

- V )2 



T J 

V t + — — 
0 

0 

OJ 

^0 

(Uo) 

L2 -( 

(V 

- 




V- - 

0 

2 A 

^0 

ihi) 


where Li is the minimum, and L2 the maximum distance to fly in interval t 
when speed is between Vq and Vf. The relation between the required speed 
profile and the distances Li and L2 is shown in figure IT 5 where is the 

aircraft speed in the constant-speed segment. The distances Lx and L2 are 
represented by the shaded areas of the speed profiles in the left-hand column 
of figure IT- The profiles to be used are in the right-hand column. The 
required profiles are selected by comparing the shaded areas of left-hand 
column profiles with those of the right. If L < Lx, the required profile 
consists of deceleration, constant speed, and acceleration; if L > L2 , the 
required profile consists of acceleration, constant speed, and deceleration. 
For Lx < L < L2, the required profile consists of deceleration, constant 
speed, and deceleration segments if > V^:*, or of acceleration, constant 

speed, and acceleration segments if Vq < V^. 

The "profile constants” are the second-segment speed Vn, starting time 
tx and ending time t2- With permissible acceleration (or deceleration) of 
ax in the first segment and as in the third segment, is computed by 

equating L with the area under the speed-history curve: 


(V •- V ^)2 (V - V .)2 

L = V t - 4 - ^ ^ 

n 2ax 2as 


(U2) 


Equation (^2) is quadratic in and can be expressed as 


Vf 

'n^2a3 ^ 2ai / ai ” a.^) \ 2as 




V 2 

o_ 

2ax 


- L = 0 


(U 3 ) 


After solving for Vj^ , the remaining profile constants tx and t2 are 


ti = 


V - V 
n o 

ax 


to = t - 


- V 
f n 

Q'3 


( 4 U) 

(h 5 ) 


The profile is physically meaningful only if tx < t2; if tx > t2 , the 
final time must be changed. The profile degenerates to two segments if 
^n “ ^o ^n “ segment if = V^. 
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The final step is to compare Vj^ with Vjo±n and V-mgy . If 
Ymin ^n Ymax* profile is valid. If Vn > Vmg y » the time permitted to 
the final point must be increased. If Vn < V-m-1* r) > the path length must he 
increased. 


Path Stretching 


This section considers a technique for increasing the path length. Path 
stretching is applied only to the straight-line portions of the tra^jectory. 

Any technique for path stretching must provide for continuous increases in 
path length, so that the difference between the length of the stretched path 
and that of the unstretched path can range from zero to any desirable number. 
The proposed technique provides for this type of continuity when the distance 
between the initial and final point is greater than four turning radii. When 
this distance is less than four turning radii, a continuous range of stretched 
path length is not always possible. 

The patterns for the stretched path consist of circular arcs and line 
segments. Although different radii may be used for different arcs without 
introducing any complexity into the path-stretching algorithm, the radii of 
all arcs in a pattern are assumed to be the same to simplify presentation of 
the technique. 

Path-stretching patterns are generated by three circles connected by tan- 
gent lines as shown in figure l8(a), to form a path confluent with the 
unstretched path at the initial and final points. Two of the circles, Ci and 
C2, are fixed and tangent to the unstretched path at Pq and P^, respectively. 
The third circle, C3, moves to generate different path lengths. The center of 
the third circle is restricted to move along a curve consisting first of a 90® 
arc with a two-unit radius and then of a half line, as shown in figure 19. 

The pattern is completed by connecting the three circles with the appropriate 
tangent lines. 


If the center of C3 lies on the half-line portion of the curve, the pat- 
tern consists of five segments and is "nondegenerate.” Otherwise, the pattern 
contains fewer segments and is "degenerate." The nondegenerate pattern is 
used most often because it provides a large range of path length. Figure I8 
shows the four patterns generated by placing the center of C3 at points P^, 
Pp, Pq, and Pp of figure 19 - As the center of C3 moves along the locus of 
centers from Pa» the difference between the stretched and the unstretched 
path length increases monotonically from zero to any positive number. Since 
the distance between P^ and P^ is not less than four radii, circles Cj, C2, 
and C3 cannot intersect. At worst, C3 is tangent to Cj and C2; thus the 
tangents required to complete the pattern can always be constructed. 


The relationship between the direct path length (the distance between 

Pq and Pf) and the stretched-path length Lg is plotted in figure 20 . The 
location of the center of circle C3 on the locus of centers (fig. 19) is a 


parameter. Between and Pq , the parameter is 0 g ; beyond 

eter is the perpendicular distance Y between the center of 

^3 


Pq , the param- 
C3 and the line 
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connecting and C2 in units of turning radii. For any fixed value of 
in the interval of [U, <») , the value of Lg increases continuously as C3 
moves along the locus of centers starting from point P^. Moreover, it can be 
concluded from the characteristics of these curves that for any given combina- 
tion of values of (L^, Lg ) , where ^ [U, «>) and Lg 2 a unique curve 

intersects point Lgl. The parameter value that generates that intersect- 

ing curve defines the solution of the path-stretching nroblem. The equation 
for the unknown value of the parameter is nonlinear in the parameter and can 
be solved by known techniques; it is derived in the appendix. 


Algorithm for Controlled Time of Arrival 


The algorithm for synthesizing controlled time -of- arrival profile con- 
sists of speed-control, path-stretching, and decision logic (fig. 2l) . The 
profile is synthesized in three steps: 

(1) Determine l^mln L^ax» ^^e minimum and maximum distances the air- 
craft is capable of traversing, for the given Vq, Vf, Ag^, A^, Vj^in 9 and 

t . 

(2) Determine if L is within the range of [Lniin^ Ljj^ax^ ’ 

L < stretch the path length to Lg so that Lg < 

proceed to step 3, (b) if L > return message stating that the given 

time t is too short for the aircraft to traverse the given distance L, (c) 
if Ljj^in - L < Ljj^ax’ Proceed to step 3. 

(3) Determine the profile pattern and compute profile constants t^, 

and t2- 


The distances I^^ax ^-^±n correspond to the areas under speed-profile 
curves whose constant-speed segments are flown at ^min^ and are com- 
puted by substituting the speed and acceleration limits into equation (U2): 


L = V t 
max max 


= V t 
max 


L . = V . t 

min min 


= V . t 
min 


(V - V )2 (V - 
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When t is too short for the aircraft to accelerate to 
corresponding maximum-distance profile consists of tvo segments : an accelera- 

tion segment from Vq to V■^p ( , and a deceleration segment from V^p to 
Yf, The upper speedy ^-upj is computed by equating the sum of the time dura- 
tion of the two segments to t: 


(t - ti) = 


V - \ ^ Vup - Vf 


= t 


Upon rearranging terms ^ 


V 

up 




( 48 ) 


The area under the curve is 


V - V 

o / 

L = V ti + ti + V„(t 

max o ^ 2 f 


V - V 

ti) + - ^ (t - ti) 


''up ''o ' 


2A 


^ 

2|a. 


( 49 ) 

(50) 


Similarly, the minimum- distance profile for small t consists of two 
segments: a deceleration segment from to (a speed greater than V-m-i n ) > 

and an acceleration segment from to V^. For 

V _ V V - V 

ti + (t - 1 1) = + 7 = t 


Upon rearranging terms , 


= 


A.J a 


( 51 ) 


The area under the curve is 


V - V V - V 

V! 

2 A, 2A 

d a 


(52) 

( 53 ) 
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Second, a comparison is performed to ascertain that the controlled time of 
arrival can he accomplished by speed control alone. If not, then the given 
nominal path length L is increased to Lg by path stretching as specified 
by k: 


L 

s 


kL + (1 - k)L . 
max min 


(5lt) 


where 0 < k < 1. Third, the path length L (or Lg) is compared with Li and 
L 2 to determine the pattera for the speed profile, and the trajectory con- 
stants are computed by equations (Us), and (^+5)* 


Time Reference 

In the controlled-time-of-arrival problem, time-to-go is an important 
variable to pilot and air-traffic controller. As arrival time approaches, air 
traffic control must know by how much this time can still be changed. The 
earliest arrival time latest arrival time achievable by 

speed control alone are 


t . 

min 


L + 


(v 

max 
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V 

max 
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max 


L - 


(V --V . )2 
o min 

2]aJ 


(V„ - V . 
f min 

2A 


min 


(56) 


These times could be used as decision variables instead of L^j^^x ^in 
the algorithm. If t < the trajectory is beyond the aircraft capabil- 
ity. If tj^±rx S t < tjj^ax 9 controlled time of arrival can be accomplished 

by speed control alone, but if t > tmax* PO-th stretching is also required. 


Vertical Guidance 

Criteria for generating the vertical profile were described earlier. The 
vertical profile is defined by the time to begin descent to the desired final 
altitude. This time is calculated by subtracting from the specified arrival 
time, t, the time intervals required to change the approach speed to the spec- 
ified final speed in level flight, and to change the cruise altitude to the 
final altitude. The first time interval is profile constant t 2 j available 
from the already calculated speed-time history, and the second interval is 
given by (h^ - hf)/ii. 
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Summary 


This section has developed procedures for synthesizing speed profiles to 
achieve specified arrival times at specified points on the approach path, and 
for modifying the length of the approach path if the speed profile required 
for a specified arrival time is outside the aircraft’s performance envelope. 

A simple technique for vertical guidance was also described. The simplicity 
of the technique resulted from the use of parameters already calculated in the 
synthesis of horizontal and airspeed profiles. The synthesis algorithms dis- 
cussed in this and the previous sections are illustrated in the next section 
by an example. 


OVERALL FLIGHT PROFILE SYNTHESIS - AN EXAMPLE 


A numerical example will be used to illustrate how a complete, three- 
dimensional trajectory with time constraints may be synthesized by applying 
the techniques developed for horizontal guidance, controlled time of arrival, 
and vertical guidance- It is assumed that an aircraft has reached the last 
waypoint before the outer marker. Wind velocity is assumed to be zero. If it 
is not, the trajectory to be synthesized must be modified somewhat to account 
for the wind velocity. 


Problem Description 

Consider an aircraft whose current position is 21.8 km (13-56 miles) from 
the OM at an azimuth of 292®, It has a heading of 2l6® , an altitude of 
1520 m (5000 ft), and a speed of l49.6 m/sec (290 knots). The flight profile 
to be synthesized must guide the aircraft to arrive at the OM six minutes 
later, aligned with the runway at a final altitude of m ( 15 OO ft) and a 
final speed of 67 m/sec (l30 knots). The acceleration and deceleration capa- 
bility of the aircraft is O. 6 I m/sec^ (2 ft/sec^), and the maximum and minimum 
speeds of the aircraft are 15^-5 m/sec (300 knots) and 67 m/sec (l30 knots). 
The sink rate (rate of change of altitude) is 305 m/min (lOOO ft/min). The 
aircraft must maintain cruise altitude (1520 m in this case) as long as possi- 
ble, descend at constant speed, and hold altitude during speed changes. A 
minimum turning radius of 6 .U 5 km (4 miles) is used during horizontal maneu- 
vers. The conditions and constraints for this example are summarized in 
figure 22. 

As outlined earlier, the complete profile of the aircraft’s time-position 
relationship is synthesized in three sequential steps: horizontal guidance, 

controlled time of arrival, and vertical guidance. Each subproblem and its 
methods of solution are discussed separately. After the complete profile is 
solved, it is summarized as a command sequence, which gives the order, initia- 
tion times, and termination times of commands required to fly the profile. 
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Horizontal Guidance 


From the description of the example, the horizontal guidance problem is 
classed as problem (a), and steps for synthesis are: (l) compute pattern 

parameters 5 (2) determine the flight pattern, and (3) compute turning angles, 
linear path length, and total path length. 


First, the pattern parameters d^ and d^ are computed using equations 
(2) and (3): 


d 


o 


yp 

R sin ijj 


5.32 


d 


f 



R tan ij; 


-5.6 


The heading of 2l6° , converted to modulus l80^ , is -l44°. Since it is nega- 
tive, a reflection about the d^ axis is required for pattern determination: 
^r = 

Second, the flight pattern is determined by comparing d^ and d^ with 
the appropriate switching points in the switching diagram (fig. 8(c);. Com- 
paring d^ with where the value of d^+ is computed by equations listed 

in table 3: 


3 1 

d + = ^ — + — = 3.73 

oa . ,r r 

sin i|j tan \\) 


so d^ > d +. The value of d 
o oa 


is then compared with d^+ where 


= 1 ^ 

fy tan ip 


sin ip 


- 2. 43 


so that d^ < These inequalities establish that the flight pattern is 

RSR and that the reflected flight pattern is therefore LSL. 


Third, the turning angles, curved path length, linear path length, and 
total path length are computed by the equations for the LSL pattern in table 

4: 


PiPo = /(x + R sin + (R - R cos + y )^ = 9*65 km (6 miles) 

^ o o 


ip = tan' 


-1 


-R + R cos ^ - y^ 

-X - R sin }p 
. o ^ , 


= 118^ 




ij; - + 360° = 98° 
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^ ^ = 11 km (6.84 miles) 


:= ijj^ = 118 ° 

^f "" \so^^ "" ^ ^ miles) 


L = P^Pi + P 1 P 2 + PaP^ = 34 km (21.1 miles) 
This concludes the computation of the horizontal trajectory. 


Controlled Time of Arriyal 


The numerical inputs for this problem are those specified in figure 22, 
with the addition of the path lengthy 34 Ion (21.1 miles), obtained in the 
preceding section. 


The logical path traced in synthesizing the speed profile is shown with 
heavy lines in figure 23. At node A, the time interval, Ata, for the maximum 
distance profile is computed to be 208.1 seconds. Since Ata is positive, the 
maximum distance profile consists of three segments. Next, Lmaxj distance 

corresponding to this profile, is computed before node ^ : 


L = V t 
max max 


(V 

max 



2A 
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= 49.1 km ( 30.5 miles) 


(v 

max 
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Since the given distance L for this problem (34 km) is less than I>niax 9 
profile synthesized is within the aircraft’s capability. 


Similarly, at node ^ , Ata for the minimum distance profile is computed 
to be 225 seconds. The positive value for Ata implies that the minimum dis- 
tance profile contains a constant speed segment. Before node ^ , the dis- 
tance corresponding to this profile is computed and then compared with L: 


L . 

min 


(V - V . )^ (V. - V . )^ 

, o min f min 

— V , "G I ■ ~ > 


min 


2 A. 


2A 


= 29.6 km ( 18.4 miles) 


Since L is greater than ^min » stretching is not necessary. 

Two further comparison tests establish a three-segment profile pattern of 
deceleration/constant speed/deceleration. At node ^ , the speed for the 


32 


constant-speed segment is computed by 'solving for in equation (Us): 


V = 86 m/sec (i 6T knots) 
n 


Finally, the values of ti and t 2 j the time instants corresponding to a change 
from deceleration to constant speed and from constant speed hack to decelera- 
tion, are computed hy substituting numerical values into equations (H4) and 

(^ 5 ): 


*^1 



103.6 sec 


t2 = t 



360 - 31. U = 328.6 sec 


This concludes the synthesis of the speed profile. 


Vertical Guidance 


The time duration, At;j^, required to descend from a cruising altitude of 
1520 m (5000 ft) to a final altitude of U56 m (1500 ft) at a sink rate of 
305 m/min (lOOO ft/min) is 




1520 - U56 
305 


= 3.5 inin = 210 sec 


Since the aircraft is permitted to descend only when the speed is constant, it 
must descend within the time interval 

[ti, tz] = [ 103 . 6 , 328.6] 

Since this time interval ( 226.3 sec) is greater than At^ , the aircraft can 
descend at constant speed. Since the aircraft is also required to maintain 
cruise altitude as long as possible, the altitude change vill take place dur- 
ing the last portion of the constant-speed time interval, from 
ti = t 2 “ At^ = 118.6 sec to t 2 - 


Overall Command Profile 


The overall command profile expresses the horizontal guidance, controlled 
time of arrival, and altitude profile in a time sequence. The time for ter- 
mination of the first left turn, t^^ , determined by the time required to 

traverse the arc length 


P Pi: 
o ^ 


P = Tpi (V 
o o 

o 


H- A^t)dt = V^tpj + 


2 ^Pl 
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(An additional equation is required if is zero over a segment of the 

turn.) Substituting the kno-wn numerical values into this equation gives the 
value of 90*53 seconds for tp^ . 

The time for initiation of the second left turn, tp 2 , is computed next. 
The distance traversed vhile the aircraft is decelerating from to Vf, ^ 52 » 
is the area under the last segment of the speed profile curve: 


Substituting the known constants gives a value of 2,k km (1.49 miles) for 
Comparing the numerical value of with the arc length 


= 13.25 km (8.24 miles) shows that t^^ occurs in the arc Then 

t 


P2 


is computed by subtracting the time to fly the constant-speed segment of 


the arc 




2 ^^ from t 2 : 


- ‘3 2 


tp2 = *2 - 


n 


= 202.23 seconds 

A complete command profile can now be constructed from a time history of the 
three profiles. Figure 24 is a complete profile of the overall trajectory. 
Figure 25 is the corresponding horizontal trajectory with points Pq , ? 2 , 

P3 9 P49 P59 Pf marked to indicate the beginning or termination of com- 
mands and the times at which they occur. Table 7 summarizes the entire com- 
mand sequence. If a ground-based computer performs the trajectory synthesis, 
table 7 gives the essential information the on-board guidance system must 
receive to fly the trajectory. 


CONCLUSIONS 


The technique developed in this report for automating terminal-area 
guidance, when properly implemented, enables an aircraft to fly from a given 
initial position, altitude, heading, and airspeed to a specified final posi- 
tion, altitude, heading, and airspeed in a prescribed time. The technique is 
based on the formulation of a class of principal subproblems in terms of which 
any flight profile can be expressed from takeoff to landing. The required 
subproblems are established by a study of air-traffic control procedures, 
instrument- flight maneuvers, and terminal-area route structures. They are: 

( 1 ) Horizontal guidance: (a) Guiding an aircraft from any initial posi- 

tion and heading to a specified final position and heading, (b) guiding an 
aircraft from any initial position and heading onto a straight line of speci- 
fied location and heading, (c) guiding an aircraft from any initial position 



and heading to a specified final position, and (d) guiding an aircraft 
smoothly from one straight line onto another intersecting straight line. 

(2) Controlled time of arrival: Guiding an aircraft (along the pre- 

selected horizontal trajectory) from an initial airspeed to a specified final 
speed in a specified time interval. 

(3) Vertical guidance: Guiding an aircraft (along the preselected hori- 

zontal trajectory and speed profile) from an initial altitude to a specified 
final altitude. 

The criterion used for selection of the horizontal trajectory is to mini- 
mize the total path length subject to the constraint imposed by a minimum 
turning radius. The horizontal trajectory consists of arcs of circles and 
segments of straight lines which form a confluent path connecting the initial 
and final points. It is synthesized in two steps: First, a flight pattern is 

assigned, then arc lengths and straight-line distances are computed to specify 
the horizontal trajectory completely. 

For the speed-profile synthesis, the horizontal distance determined 
already is compared with the maximum and minimum distances the aircraft can 
traverse (for the specified initial speed, final speed, time of arrival, maxi- 
mum and minimum allowable speed, and acceleration and deceleration capabilitv). 
If the horizontal path length is greater than the maximum distance, no speed 
profile is possible. If it is between the minimum and maximum distance, speed 
variation alone can control time of arrival. If it is below the minimum dis- 
tance, path stretching of the horizontal trajectory is necessary, and the con- 
trolled time of arrival is achieved by a combination of speed control and path 
stretching. 

When the horizontal trajectory and the speed profile are determined, the 
vertical guidance trajectory is synthesized. The criterion for this synthesis 
is maintaining cruise altitude as long as possible consistent with the speci- 
fied descent rate and with the required final altitude. 

A composite command profile is generated by superimposing the already 
determined horizontal, speed, and vertical profiles. Essentially, this com- 
mand profile is a time history of commands (such as decelerate, turn left, and 
descend) to be issued at the proper time. This profile can provide automatic 
guidance of an aircraft for transitioning from any initial state to any speci- 
fied final state, where a state is specified by the aircraft coordinates, 
altitude, speed, heading, and the clock time. 

In general, the overall flight profile of an aircraft, starting from 
takeoff to landing, may be decomposed into a sequence of waypoints. The 
flight profile for each pair of adjacent waypoints can readily be synthesized 
by applying the various guidance problems discussed. Therefore, the overall 
flight profile, starting from takeoff to landing, may be synthesized by con- 
catenation of flight profiles for consecutive pairs of waypoints. 
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The advantage of this technique lies in its versatility and efficiency. 

It requires no precomputation of trajectories, and can provide fast, in-flight 
changes in trajectory as required. Since the only aircraft parameters used in 
synthesizing the trajectories are performance limitations, the technique is 
applicable to any aircraft. Finally, by synthesizing trajectories subject 
only to aircraft performance limitations and ATC constraints , the technique 
uses airspace efficiently and therefore makes possible reduction of separation 
requirements and thus achieves higher landing rates . 

Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 9^035, February 11, 1972 
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APPENDIX 


COMPUTATION OF PATH LENGTH FOR STRETCHED PATH 


For the coordinate system used to compute the stretched-path length, the 
origin is^at the center of circle Ci and the positive x-axis is in the direc- 
tion of Hq; the unit of distance is the turning radius. For the nondegener- 
ate pattern in figure 26, the parameter is and the total length of the 


stretched path is 


^C3 


^s " ^0^1 ^1^2 + P2P3 + P3P4 + 


where 


(Al) 


P P = — 
o 1 2 


(A 2 ) 


P1P2 ^C3 


P2P3 ~ 2 ^ 


PJ5T = - ("o' - - - 


P4P^ = ^ 


(A 3 ) 

(aU) 

(A 5 ) 

(a6) 


The angle X is given hy 


X = 7T - tan 


L - 2 

-1 o -1 

- cos ^ 

>^C3 


4^3 * (L„ - 2)^ 


(AT) 


Equations (A5) and (AT) are derived hy solving for parts of right triangles 
P5C3C2 and C3P3P^, By adding expressions (A2) through (A6) we obtain the 
expression for the total path length. 


_ L -2 

L = 3 tt + y^ - 2 tan”^ — — - 2 cos"^ 
^ "3 ^C3 


^ + M + (L.-2)2 - u 

42 + (L -2)^ ^3 


C3 


(A8) 


For the degenerate pattern in figure 2T, the parameter is 6g and the total 
path length is 


+ P1P3 + P3P4 + P4P, 


(A 9 ) 


where 


PqPi = 0s 


(AlO) 
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(All) 


P1P3 = 63 + A 

P3P4 = A cos^ 0 + (L - 2 sin 9 )^ - ij- 

^ ^ SO s 


= X 

Solving for parts of right triangles C 3 P 5 C 2 and C 3 P 3 P^, 


X = tan 


L - 2 sin 0 

-TO s -1 

^ a - COS 

2 COS 0 


/(L - 2 sin 0 cos 0^^ 

a c 


The total path length can then he written in terms of X hy adding (AlO) 
through (AI 3 ): 


L =20 + 2X + A cos 2 e + (L - 2 sin 0 )^ - 4 

s s so s 


(A12) 

(A13) 


(Alh) 

(A15) 
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TABLE 1.- SUMMARY OR TERMINAL-AREA GUIDANCE PROBLEMS 


Class 1 Brief title 

Description 

Onerational constraints 

Preferred 

synthesis criterion 

Horizontal 

guidance 

Problem (a) 

Fly from any initial nosition 
& heading to a specified 
final position & heading 

Maxim-um bank angle 

Minimum nath length 
to achieve final 
position & heading 

Problem (b) 

Fly from any initial nosition 
& heading onto a line of 
specified location & direction 

Maximum bank angle 

Minimum path length 
to achieve final 
line 

Problem (c) 

Fly from any initial position 
& heading to a specified final 
position vith arbitrary final 
heading 

Maximum bank angle 

Minimum path length 
to achieve final 
point 

Problem (d) 

Transition between intersect- 
ing straight lines 

Maximum bank angle 

Smooth transition, 
minimum overshoot 

Path 

stretching 

Increase length of flight 
path by a specified amount 

Maximum bank angle 

Unknown 

Speed 

control 

Controlled 
time of 
arrival 

Determine speed profile to 
achieve specified arrival 
time at a specified location 
with specified final speed 

Minimum & maximum 
sneed & maximum accel- 
eration & deceleration 

Simplicity of speed 
profile 

Vertical 

guidance 

Descent 
from cruise 

Determine flight path in ver- 
tical plane to achieve a 
1 specified final altitude at a 
: specified position 

Maximum descent rate 

Maintain cruise 
altitude as long as 
possible & descend 
at constant h 




TABLE 2.- OPTIMAL FLIGHT PATTERNS FOR MINIMUM-ARC-LENGTH 

HORIZONTAL GUIDANCE 



Additional patterns required 
for optimum solution (ref. 9) 
(not used) 

Maneuvers* Geometric shape 




* Significance of symbols in column 
as follows: 

L — left turn 
S — straight flight 
R — right turn 
Pq ~ initial position 
Pf — final position 










TABLE 4.- TRAJECTORY PARAMETER FORMULAS FOR PROBLEM (a) 


j 1 't's ’ I *^1 ’ ; ’ 1 

; Pattern; ?lP2» Heading of straight-line Heading change i Heading change > 

' type “ Length of straight-line segment segment > in initial turn. in final turn i 


RSR 

;/(x-R sini|^)^+(R-R cosijj-y)^ 

R-R cosi^-v 
arctaxL . , 

-x+R sini|/ 

>.=l 

Eh 

1 ! 

2ir-i|)g :'|<g>0 

0 ;1.3=0 ; 

RSL 

2R___ „hprr 

'>3=. 

^1+^^ ;\(^l+ip^ <2 tt 

1 

1 


( 

'/'f='l'g 

, ■ ■, , wnere 

tan 

i|jl+i|;' -2 tt ;\);i+il;’^2'fT 



■ ■ ^ ■ - ■ - - ■ - ^ , , -R-R cosil^-y l 

/(-x+RsiniJj)^+(R+Rcosi|)+y)^ ^ ^ -x+R sintj; j| 


LSL 

/(x+R sinif^)^+(R-R cosi)j+y)^ 

j 

-R+R cosij^-y 

arctan 5 — : — 1 ^ 

-x-R siwp 

%=■ 

_ip-i|;3+2ii 

4'f='l's 

LSR 

2R 

. , , where 

tan 

... 2R 

i|j*=arcsin 

/( x+RsiniJj ) (R+Rcos i|)-y ) ^ 

^s=. 

whe] 

iJjj-ij)’ ; >0 

i|)l-t|)+2iT <0 

R+R costj^-y 

re ibi=arctan _ . , 

-x-R sini^ 

4»i=( 

'I'-'I's ; fe’I's 

tp-)|;^+2'Tr;i|;<i|;g 

2fr-tpg;tjjg>0 

'Pf= 

0 ; <J's=0 






TABLE 5 .- TRAJECTORY CONSTANTS FOR PROBLEM (b) 


Pattern 



Range of 

type 

i Specification of region 

Trajectory constants 

turning angles 


-7T < < -tt/2 

d ^ 

, \p^ = 3tt/ 2 + ij; 
ipf = tt/2 

t = -d - 1 + cos 4’ 

0 ^ ^ TT 

ir/2 S ip ^ TT 

\pj_ = ip - tt/2 

. ijjf = tt/2 

d ^ S3U) 

C\J 

II 

1 



il = -d - 1 + cos ip 



-tt/ 2 ^ ff/2 

d ^ S2(ii^) 


tt/2 < ^\) <tt;S2(i|^) < d < Si(\J;) 


0 < ij; < tt;S 3 (i|;) < d < S 2 (^) 
or 

-tt/2 < }\) < 0 ;Si(\p) < d < S 2 (^) 


-if 1 + cos ^ - 


T|;i = -cos" 


, _ 1 / 1 + cos \j; 

= -cos 2 “-^ 


, -1 / l + cos 1 |^ - d\ . , 

\l)l = cos j + \\) 


\\)^ - cos 


-iILjL 


cos ij; - d 


+ ^ , 0 < < tt/2 

+ 2tt ’ 7T < < 3tt/2 

0 ^ ^ 3 V 2 

, 0 < < tt/2 










TABLE 6.- TRAJECTORY CONSTANTS FOR PROBLEM (c) 



1+5 



TABLE 7.- OVERALL COMMAOT) PROFILE 


Time , 

Point second 


Polar coordinates* 


Commands 


0 21.8 (13.56)5 292 Begin left turn and 

begin deceleration 

90.5 18 (11.18), 266 Fly straight 

103.6 17.2 (10.58), 263 Hold speed 

118.6 15.95 (9.9), 261 Begin descent 

202.2 11.3 (6.85), 239 Begin left turn 

328.6 2.38 (I.U8), 191 Hold altitude and 


191 Hold altitude and 
begin deceleration 

0 Fly straight and 
hold speed 


*Distances in km (miles), angles in degrees 





I 

I 

I Oakland 



Altamont*^ 

/ 

/ 

/ 






\ / 

^-/ 

Woodside VO^ 


\ 


\ 

\ 

\ 



Area of variable 
route structure 






,^San Jose VOR 


n mi. 


— I 



Figure 1.- SPO terminal-area approach routes. 
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Figiore 6 


(a) between P - and P - ; Pf between P,± and P,- and above P,- 

o oa oa T fo f7 fa 



(c) P^ between P - and P + ; Pf between P,+ and Pfc 
o oa oa ' iy 'o 


- Switching points and near distance trajectories if 
and P-p lies between P_. and P 4-- 


p + 


lies above 











Figure 9*~ Switching 




■A — A — A — — ^A — A — A — ^A — ^A — A 



(b) l// = 7T 

uns for ij; = 0° and ij; = ir radians. 























a 

CO 


Figtire l4.- Trajectory algorithm for problem (b). 
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Figure 19*" Locus of centers for c 


Locus of C3 



Pf nf 



3 * 
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Figure 20.- Unst retched versus stretched path length. 
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Conditions 

Initial 

Final 

Constraints 

X coordinate of Pq, km 

5.2 

0 

Turning radius 

6.45 km 

y coordinate of Pq, km 

-12.5 

0 

Acc., dec. 

0.61 m/sec2 

Aircraft heading, degrees 

216 

0 

Max. speed 

154.5 m/sec 

Speed, m/sec 

149.6 

67 

Min. speed 

67 m/sec 

Altitude, m 

1520 

456 

Sink rate 

305 m/m in 

Time, sec 

0 

360 




Figiire 22.- Parameters for exaaple problem. 
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